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Abstract: 

We carry out a detailed calculation of the cross section of pion diffraction dissociation into 
two jets with large transverse momenta, originating from a hard gluon exchange between 
the pion constituents. Both the quark and the gluon contribution are considered and in 
the latter case we present calculations both in covariant and in axial gauges. We find that 
the standard collinear factorization does not hold in this reaction. The structure of non- 
factorizable contributions is discussed and the results are compared with the experimental 
data. Our conclusion is that the existing theoretical uncertainties do not allow, for the 
time being, for a quantitative extraction of the pion distribution amplitude. 
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1 Introduction 



After thirty years of Quantum Chromodynamics, many aspects of hadron structure re- 
main poorly understood. The bulk of the existing experimental information comes from 
parton distributions that can be interpreted as one-particle probabilities to find quarks, 
antiquarks and gluons carrying certain momentum fractions of the parent hadron. The 
parton distributions are inclusive quantities, in the sense that contributions of different 
parton states are summed over, and therefore provide us only with global information 
about the hadron wave function. It is well known that hard exclusive reactions allow, 
in principle, to separate contributions of different Fock states and study the momentum 
fraction distributions of the components with the minimum number of Fock constituents 
at small transverse separations, dubbed distribution amplitudes. In practice, however, 
progress in this direction had been limited, due to both experimental difficulties to isolate 
exclusive amplitudes and theoretical problems to formulate their quantitative description. 

Classical applications of QCD to exclusive reactions addressed electromagnetic form 
factors JT], 0, [!|. Nowadays it is almost generally accepted, however, that the asymptotic 
behavior of form factors is only achieved for very large momentum transfers. In addi- 
tion, the form factors involve a convolution (overlap integral) of distribution amplitudes 
which makes it difficult to extract them directly from the data. In recent years, the list of 
applications of QCD factorization to hard exclusive reactions have been increased signifi- 
cantly, with hard exclusive meson production || |4| and deeply-virtual Compton scattering 
(DVCS) |5], H providing notable examples. This raises hopes that suitable hard pro- 
cesses can be found in which hadron distribution amplitudes can be studied in a more 
direct way. 

In particular, coherent diffract ive production of dijets by incident pions (or photons) 
on nuclei has attracted a lot of attention. The subject of diffraction is very old. Ever 
since the classic works in early 50' on the diffraction breakup of deuterons |J it was 
known that the momentum distribution of the proton and the neutron in the final state 
is close to their momentum distribution as deuteron constituents. More recently, the 
same idea gave rise to the method of so-called Coulomb Explosion Imaging |IJ which 
is widely used to study the high-momentum tails of wave functions of small molecules. 
To our knowledge, the pion (and photon) diffraction dissociation into a pair of jets with 
large transverse momentum on a nucleon target was first discussed in [fTI |. In Ref. Jl2|] 



the possibility to use this process to probe the nuclear filtering of pion components with 
a small transverse size was suggested. The A-dependence and the q\- dependence of the 



coherent dijet cross section was first calculated in |13 |. In the same work it was argued that 



the jet distribution with respect to the longitudinal momentum fraction has to follow the 
quark momentum distribution in the pion and hence provides a direct measurement of the 
pion distribution amplitude. Recent experimental data by the E791 collaboration |TJ], [Tj| 
indeed confirm the strong A-dependence which is a signature for color transparency, and 
are consistent with the predicted ~ 1/g^ dependence on the jet transverse momentum. 
Moreover, the jet longitudinal momentum fraction distribution turns out to be consistent 
with the ~ z 2 {l — z) 2 shape corresponding to the asymptotic pion distribution amplitude 
which is also supported by an independent measurement of the pion transition form factor 
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Figure 1: Kinematics of the coherent hard dijet production ir — ► 2 jets. The hard scattering 
amplitude Th contains at least one hard gluon exchange. 



7777* [16|] 



After these first successes, one naturally asks whether the QCD description of coherent 



dijet production can be made fully quantitative. Two recent papers [|17] and [[Tj| address 
this question, with contradictory conclusions. In the present work we attempt to clarify the 
situation and develop a perturbative QCD framework for the description of coherent dijet 
production that would be in line with other known applications of the QCD factorization 
techniques. Some of the results reported in this paper were published earlier in a letter 



format [T9|]. An approach close to ours was suggested independently in pOfl , |2T| . We 
disagree with |]20| , |2l| on several important issues and the reasons for this disagreement 
will be elucidated in what follows. 

The kinematics of the process is shown in Fig. |I[ For defmiteness, we consider 7r~ 
scattering from the proton target. The momenta of the incoming pion, incoming nucleon 
and the outgoing nucleon are pi,p2 and p' 2 , respectively. The pion and the nucleon masses 
are both neglected, pi = 0, p\ = (p'2) 2 = 0. We denote the momenta of the outgoing quark 
and antiquark (jets) as qi and g 2 , respectively. They are on the mass shell, qf = q\ = 0. 
We will use the Sudakov decomposition of 4-vectors with respect to the momenta of the 
incoming particles p\ and p 2 - For instance, the jet momenta are decomposed according to 

rp 1 rp 

qi = zp x + — p 2 + qi± , q 2 = ZP1 + — P2 + <?2± (1-1) 

zs zs 

such that z is the longitudinal momentum fraction of the quark jet in the lab frame. We 
will often use the shorthand notation: u = (1 — u) for any longitudinal momentum fraction 
u. The Dirac spinors for the quark and the antiquark are denoted by u(qi) and v (q^)- 

We are interested in the forward limit, when the transferred momentum t = {P2—P2) 2 
is equal to zerd^ and the transverse momenta of jets compensate each other qu_ = qj_, 
q 2 ± = —q±- In this kinematics the invariant mass of the produced qq pair is equal to 
M 2 = q\j zz. The invariant cm. energy s = (p\ +P2) 2 = 1p\P2 is taken to be much larger 

Hf the target mass m is taken into account, the momentum transfer t = (j>2 — p'2) 2 contains 
a non- vanishing longitudinal contribution and is constrained from below by \t\ > to, where to = 
(m 2 M ) / (s — m 2 ) 2 , M 2 being the invariant mass of the dijet. 



3 



than the transverse jet momentum q±. In what follows we often neglect contributions to 
the amplitude that are suppressed by powers of q\/s. 

From the theoretical point of view, the principal question is whether the relevant 
transverse size of the pion r± remains small, of the order of the inverse transverse momenta 
of the jets r± ~ 1/<?_l- In this paper we investigate the possibility that the amplitude of 
hard dijet coherent production can be calculated using the standard collinear factorization 
in the form suggested by Fig. [I]: 

1 1 

Mn-^jets = ^ dz' dx^^Z^/J^Tfjiz^Xi^p) F^Xi^/Jp) . (1.2) 

p=q,q,g o o 

Here ^(z 1 , up) is the pion distribution amplitude, and F^(xi,fip) is the generalized 
(skewed) parton distribution p = q,q,g [0, |23|, in the target nucleon or nucleus; 
X\ and X2 = X\ — ( are the momentum fractions of the emitted and the absorbed partons, 
respectively. The asymmetry parameter ( is fixed by the process kinematics: 

( = M 2 /s = ql/zzs. (1.3) 

Tjj(z' ,Xi, up) is the hard scattering amplitude involving at least one hard gluon exchange 
and [ip is the (collinear) factorization scale. By definition, the pion distribution amplitude 
only involves small momenta, k± < fiF, and the hard scattering amplitude is calculated 
neglecting the parton transverse momenta. In this paper we present a detailed calculation 
of the leading-order contribution to T H (z', Xi, up) corresponding to a single hard gluon 
exchange. 

We consider both the quark and the gluon contribution to the amplitude, and in both 
cases find that the corresponding hard kernels Tjj, Tfj diverge as 1/z' 2 and 1/z' 2 in the 
z' — > and z' — > 1 limit, respectively. This implies that the integration of the pion 
momentum fraction diverges at the end-points and the collinear factorization is, therefore, 
broken. Physically this means that the approximation of neglecting the incident quark 
transverse momenta becomes insufficient close to the end points, similar to what happens 
e.g. in the heavy quark limit in B-decays p5 ], albeit for a different reason. One may 



argue that Sudakov-type radiative corrections should suppress the end-point contributions 
and try to develop a modified factorization framework, by the resummation of soft gluon 
emission to all orders PS| . Alternatively, the end-point singularities may be softened on a 
nuclear target because of the color filtering of configurations with a large transverse size. 
Both possibilities are interesting and require detailed studies that go far beyond the tasks 
of this paper. We will rather assume that the end-point behavior can be tamed in some 
reasonable way, and examine the consequences. 

In particular, we will analyze the structure of the hard gluon exchange contribution in 
some detail. We will find that the structure of this exchange is such that it generates an 
enhancement by a logarithm of the energy in the region \z' — z\ <C 1. If only this logarith- 
mic contribution is retained, the collinear factorization is restored and the longitudinal 
momentum distribution of the jets to this accuracy indeed follows the shape of the pion 
distribution amplitude The hard gluon exchange can in this case be considered as 
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a part of the unintegrated gluon distribution, as advocated in [|HJ]. Beyond the leading 
logarithms in energy this proportionality does not hold. Remarkably enough, we find that 
the longitudinal momentum fraction distribution of the jets for the non-factorizable con- 
tribution is calculable, and turns out to be the same as for the factorizable contribution 
with the asymptotic pion distribution amplitude. We also find, in agreement with |21 



that the quark contribution is significant in the energy range of the E-791 experiment and 
present a new data analysis including all contributions. On the technical side, we present a 
detailed study of the light-cone limit of the relevant amplitudes, the structure of different 
absorption parts and the gauge- dependence. We believe that some of this discussion is of 
general interest and relevant for all exclusive processes in the diffractive kinematics. 

The presentation is organized as follows. In Sec. 2 we consider the quark contribution 
to the coherent dijet production, which is simpler than the gluon contribution as it does 
not involve subtleties related to the choice of the gauge. We find an end-point singularity in 
the integration over the pion quark (antiquark) momentum fraction and trace its origin to 
pinching of the integration contour over the partonic momentum fraction in the so-called 
Glauber region. We explain this result by analyzing the structure of different dispersion 
parts of the amplitude. The gluon contribution is considered in Sec. 3. We point out that 
there is a subtlety in the definition of the generalized gluon distribution, and explain how 
this subtlety is resolved in a more familiar case of hard electroproduction of vector mesons. 
Next, we present a calculation of the imaginary part of the amplitude of the coherent dijet 
production from a quark, in which the above-mentioned problem is avoided. Finally, we 
perform the calculation of the full amplitude (including the real part) using axial gauge 
and not assuming the light-cone dominance from the beginning. Similar as in the quark 
case, we observe pinching of singularities between soft gluon interactions in the initial and 
the final states. A detailed numerical analysis and the comparison to the E791 data is 
presented in Sec. 4, while Sec. 5 is reserved for the conclusions. 

One technical remark is in order. Calculations in the main text are presented using 
Radyushkin's definitions || ^] for the generalized parton distributions. The relation to 
the more commonly used symmetric notation by Ji [21 is explained in App. A and in 
App. B we collect all the final expressions in the symmetric notation. 



2 The Quark Contribution 
2.1 The Calculation 

We begin with the calculation of the quark contribution to the hard dijet production, 
which is simpler than the gluon contribution as it does not involve subtleties related to 
the choice of the gauge. For definiteness we consider 7r~-proton scattering. The quark 
contribution starts at order 0(a 2 s ) and can be decomposed in three topologically different 
contributions shown in Fig. |2], Fig. [5] and Fig. [|. In what follows we refer to them as 
the w-quark annihilation, the ci-quark exchange, and the gluon exchange (flavor-singlet) 
contribution, respectively. All three contributions are separately gauge invariant and can 
be calculated with quarks on the mass shell. In this section we work in the Feynman gauge. 
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a) 



b) 



c) 



d) 



e) 



f) 



I) 



Figure 2: The u-quark (annihilation) contribution to the coherent hard dijet production 
7r~ — > 2jets, see text. 



The antiquark contributions do not require a new calculation and can be obtained from 
the corresponding quark contributions by obvious substitutions, see below. The reaction 
kinematics and the notation for the momenta are specified in the Introduction. We define 
the generalized quark (antiquark) distribution as the matrix element of the light-cone 
operator |], |23| 



(P2W) H(y)\P2) y *=o = u(p' 2 ) MV2) ■ I d Xl [e-^^J^M - e™tov)j*( Xx 



(2.1) 



where q = u,d, . . ., etc. Hereafter we use the notation 

x 2 = x x - (, 



(2.2) 



where the asymmetry parameter ( is defined as p 2 —p' 2 = (P2 and is fixed by the kinematics 
of the reaction, see (|1.3| ). The relation to other common parametrizations of the generalized 
parton distributions is discussed in Appendix A. In turn, the pion distribution amplitude 
is defined as [§ @, § 



(o|d(o) My)K(pi))y*=o = i{viy) U / dz' e - iz '^^{z') 



dz' 1 



(2.3) 



where f w = 133 MeV is the pion decay constant. In both cases, ( |2.1| ) and (|2.3|), the 
insertion of the path-ordered gauge factor between the quark field operators is implied. 
Both distributions depend on the factorization scale fip which to the leading logarithmic 
accuracy has to be taken of the order of the transverse momentum of the jets. 

As an example, consider the first diagram in Fig. 0. The corresponding contribution 
reads: 



iM 



2d 



-1 u 



2N r . 



x- 





U{qi)l^ ^i7m(C j>2~ j2)lv j>2lvV{(l2) 

[(CP2 - q2) 2 + ie][(qi - ^Pi) 2 ][(<?2 + x 2 p 2 



\2l ' 



(2.4) 
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where u(qi) and v (q 2 ) are Dirac spinors for the outgoing quark and antiquark jets, respec- 
tively, and i, j are the jet color indices. The color state of the target quark is averaged 
over. Note that we only show the Feynman +ie prescription for one of the propagators, 
since the other two are strongly off-shell. 
A simple calculation yields 

1 1 

■ M 2a= Ar2 % - V 1 - C^(gi)75— v{q 2 )5ij dz'cp^z') dx 1 ^(x 1 )-I 2 a{z,z',x 1 ,x 2 ) (2.5) 
N c q ± s 1 1 

with the coefficient function 







2C 2 z 

h a (z,z',x 1 ,x 2 ) = — — ( F , . x , (2.6) 

Z Z\X 2 + 16) 

where Cp = (N 2 — 1)/ (2N C ), N c is the number of colors. Calculation of the other diagrams 
in Fig. is equally simple. The corresponding coefficient functions are 



I 2b (z,z',xx,x 2 ) 
I 2c {z,z',x 1 ,x 2 ) 
I 2d (z,z',xx,x 2 ) 

I 2e (z,z',X 1 ,X 2 ) 

I 2f (z,z',x 1 ,x 2 ) 



2C\z _ 
ie]' 

Cpz((z — x\) 



z'z'lxi + ie] 



N c z'z'[xi + ie][xi(z — z') — zz'( + ie] 

Cpz((z' — x\) 
N c z'z[x 2 + ie][xi(z' — z) — z'z( + ie] 

C F N c z( Xl + ((l-z-z')) 
z'z'z\x\ + ie] [x 2 + ie] 



2C F Qz 



z'lxi + ie][x 2 + ie] ' 



j {7 j r r x C F ( 2 z(l -z-z') 

l 2g [Z, Z , Xi, X 2 ) — — — — — — jr — 77 . n , \&-() 

N c z [x\ + ie] [x 2 + ie] [X\ [1 — z — z ) — zz Q + ie] 

so that the total contribution of the w-quark annihilation is given by the expression in 
Q2.5p with the coefficient function 

I{z,z\x u x 2 ) = ^ h(z,z',x 1 ,x 2 ) (2.8) 

i=2a,2b,...,2g 

instead of I 2a . 

The (i-quark exchange contributions in Fig. || can be obtained from the corresponding 
diagrams in Fig. |2| by the substitution: x\ — > — x 2 , x 2 — > —x±, z — > z, z' — > z' and an overall 
minus sign. Similarly, the antiquark (i-annihilation and ^-exchange coefficient functions 
can be obtained from the above expressions by the substitution X\ — > — x 2 , x 2 — > —x% and 
changing the overall sign. Altogether, we obtain four coefficient functions: 

I u (z, z', Xi, x 2 ) = I(z, z', x t ,x 2 ) , 
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a) b) c) d) 




e) f) g) 

Figure 3: The d-quark (exchange) contribution to the coherent hard dijet production 

7r~ — * 2jets, see text. Note that the quark and the antiquark lines entering the pion blob are 

interchanged compared to Fig. 0. 




a) b) 



Figure 4: The flavor-singlet quark contribution to the coherent hard dijet production it — ► 2 jets. 



I d {z,z' ,x u x 2 ) = -I(z,z,-x 2 ,-x 1 ), 
Iu(z,z',x 1 ,x 2 ) = -I(z,z',-x 2 ,-x 1 ) , 

Ij(z, z', x 1 ,x 2 ) = I(z, z, x 1 ,x 2 ) . (2.9) 

Last but not least, we have to take into account the diagrams in Fig. |] corresponding 
to the flavor-singlet two-gluon exchange. Note that the virtuality of both gluons is large, 
in fact equal to the transverse momentum of the jets, and, therefore, these diagrams are 
part of the quark coefficient function rather than the gluon one. A simple calculation 
yields 

C F (C(1 + z - z') - 2x x ) 



h b (z,z' , ,x li x 2 ) 
Denoting the sum of them as 



z'z'C[xi(z - z') - zz'C, + ie) 

C F ({(l-z + z')-2x 1 ) 
z'z'([xi(z' - z) - z'zC + ie] 



(2.10) 



+ hb (2.11) 

we obtain the final answer for the leading-order quark contribution to the dijet production 
in the form 

l l 

N 2 c ql 



M 



quark 



■y/l - (>(gi)75— v(q 2 )5ij / dz'(j) % (z') / dx x 



x y^-TTQpi) I&*e{z,z',Xi,X2) + I p {z,z',Xi,X2) , (2.12) 
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->-H — 

o c " i o c " i 

Figure 5: Deformation of the integration contour in the complex x\ plane in the Glauber region, 
see text. 

where the summation goes over all possible quark-parton species: p = u,u, d, d, etc. 
2.2 The End-Point Behavior 

Next, we have to examine the behavior of the coefficient functions I g i ue and I p at the end 
points z' — > and z' — > 1 and check the convergence of the integration in Q2.12Q over the 
quark momentum fraction z' in the pion. Convergence is necessary for the self consistency 
of the collinear factorization approach as it ensures that quark transverse momenta in the 
pion can be neglected. 

Consider the gluon-exchange contribution I g i ue first. In this case it is easy to see 
that Igi U e ~ V- 2 ' and I g i ue ~ 1/z' at z' — > and z' —>■ 1, respectively. Since the pion 
distribution amplitude vanishes linearly at the end-points, at least at a sufficiently high 
scale, the corresponding contribution to the result in ( [2 . 1 2j ) is finite. 

From Eqs. (|2.6| ) and ( |2T7|) it is seen that the contributions of the w-quark and d- 
antiquark are finite as wellQ. For the d-quark and the w-antiquark contributions the situ- 
ation is different, however. In this case we find 

I d (z,z',x 1 ,x 2 )\ z ^ = 2 -^^8{x 2 ) + O(i) , (2.13) 

N c zz z' 

and 

r(z } z',x 1 ,x 2 )\ z ^ 1 = ^L^5(x 2 )+0(-). (2.14) 

N c zz z' 

Note that: a) the ~ 1/z' 2 (~ l/(z') 2 ) behavior leads to the logarithmically divergent 
integral over the pion quark momentum fraction; b) this contribution is purely imaginary, 
and c) this contribution is proportional to ^(^2) = 5{x\ — (), i.e. it is due to parton 
configurations with vanishing longitudinal momentum fraction of one of the quarks in 
the target proton. In the terminology of |27| this corresponds to the Glauber region. In 
the remaining part of this Section we explain this result for the example of the d-quark 
contribution in some detail. 

Explicit expressions for the diagrams in Fig. |3| can be obtained from Eqs. ( |2.6| ), ( |2.7| ) 
using the substitution rules in (|2.9|). One can easily convince oneself that the only diagrams 
capable of producing a ~ 1/z' 2 singularity are those in Fig. ^|c and Fig. |3|g. In particular, 
the diagram in Fig. |3|c gives: 

j (7 j r T x C fK* + x 2 ) , , 

l3c{Z, Z , Xi, X 2 — — — Tf -, 77 ,_ . 1 • l^-J-^J 

N c z z [x 2 — ie\ [x 2 (z — z ) — z z( + ie\ 



■''We remind that the flavor identification refers to the particular case of n meson scattering. 
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It is seen that for finite values of x 2 this contribution is proportional to 1/z' to the first 
power. However, as z' — > the two poles in the Feynman denominators in Eq. (|2.15|) 
produce a pinch at x 2 = ( X\ = Q in the integral over x\, see Fig. [|. This pinch is 
responsible for the singular behavior of /3c, ~ 1/z' 2 . To see this, note that if the two 
poles would lie on the same side of the integration contour, the singularity at x 2 = could 
be avoided by the contour deformation, following the arguments of Ref. |28| (see also 0). 



In the presence of the pinch one can also deform the contour, but in this case an additional 
contribution arises 2 ' C 

I$° le (z, z', X!,x 2 ) = F —5(x 2 ) . (2.16) 

i v c z 

This pole contribution is entirely responsible for the leading asymptotics of in the 
z' — > region. The calculation of the diagram in Fig. §5 is very similar. In this case also 
one can deform the contour as it is shown at Fig. 5 and the leading asymptotics of I^ g is 
again given by the pole contribution 



2mC F z 2 r ^ 

'3 9 {Z,z ,x 1} x 2 ) = —- a°( x V ■ ( 2 - 17 ) 

i V e z z 



rpole 



The sum of these pole contributions gives the result in ( |2.13| ). For the w-contribution 



again, the pinch at x\ = ( appears at z J — > 1 and leads to the result (|2.14|) . 



The pinching of integration contours in the Glauber region indicates a serious problem 



with collinear factorization. It is known |29] that such pinches generally occur between 



initial and final state interactions involving soft particle (gluon) exchanges. In the case of 
the Drell-Yan production the pinches disappear in the sum of all Feynman diagrams for the 
cross section, which is due, in physical terms, to cancelation of the final state interactions 
|28|| (see also |30[). In our case the pinch occurs in the leading-order contribution and for 



quark exchange, which is unusual. Note that the problem with pinching singularities in the 
Glauber region is in principle unrelated to the end-point behavior of the pion distribution 
amplitude and is more general. We will return to this discussion in the next section. 

Our result in (|2.12| ) differs from the expression for the quark contribution obtained in 



|2l| by the above discussed pole terms. The authors of [^1|] have restored the imaginary 
part of the coefficient function by requiring that the physical amplitude only depends on 
s + ie and hence on ( — ie (with the real integration variable Xi). It is easy to see that with 
this prescription all singularities in the denominators are below the integration contour 
and do not obstruct the analytic continuation. In this sense, this prescription corresponds 
to the "true" light-cone contribution. The argument is not correct, however, since apart 
from the usual s-channel discontinuity the amplitude in question has a discontinuity in 
another invariant variable, M 2 , the invariant mass of the jets. The full imaginary part is 
only restored in their sum and the existence of two different types of cuts is reflected in the 
structure of Feynman ie prescription in the propagators. To explain this point, consider 
the different dispersion parts corresponding to the diagrams in Fig. [|c and Fig. |3]c, see 
Fig.|. 

It is obvious that the cut diagrams denoted as 2c-I and 2c-III describe the discontinuity 
in the invariant energy s, and the cut diagram 2c-II corresponds to the discontinuity in the 
M 2 -variable. Note that in the light-cone limit the s-channel cut 2c-I and the M 2 -channel 
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3c-I 3c-II 
Figure 6: Different dispersion parts corresponding to the diagrams in Fig. |2|c and Fig. ^c. 



cut 2c-II both occur at x\ = 0, see Q2.7|) . Moreover, their contributions exactly cancel 
each other. In the usual language this cancelation manifests itself in the following way. 
The part of the coefficient function, l 2c {z, z' , x±, x 2 ), which corresponds to the diagram in 
Fig. 0c has a pole at x\ — > 0. This pole does not contribute, however, to the imaginary 
part of the amplitude since the generalized quark distribution function vanishes at this 
point, J-q(x)\ x ^o — * |23], [/]]. As the result, we are left with the single s-channel cut 2c-III 
corresponding to the [xi(z-z') — zz'Q + ie] denominator in I 2c (z, z' ,x±,x 2 ). Note that the 
Feynman ie prescription to go around the pole is in agreement with the above statement 
that it corresponds to a singularity in the s-channel: The sign of ie can be understood as 
the substitution s — > s + ie, or ( — > ( — ie. 

Now let us turn over to the dispersion parts of the diagram corresponding to Fig. |3]c. 
In this case there are two possible cuts 3c-I and 3c-II that describe the singularities in the 
M 2 - and s-channels, respectively. The corresponding coefficient function I^ c {z, z', xi, x%) 
contain two poles, see eq. Q2.15Q . The [x 2 — ie] = [x\ — £ — ie] pole is related to the M 2 - 
channel discontinuity and its ie prescription corresponds to the substitution M 2 M 2 +ie, 
or £ —>■ £ + ie. On the other hand, the ie prescription of the second denominator [x 2 (z — 
z' — z')z + ie] can be understood as the substitution s — ► s + ie, or ( — > ( — ie, in agreement 
with the interpretation that this singularity occurs in the s-channel. To summarize, we see 
that for the case of the diagram in Fig. ^|c both dispersion parts contribute in a nontrivial 
way and pinching of the integration contour described above in fact occurs between the 
discontinuities in different channels. It is also seen that the pinching occurs between the 
soft quark exchange in the initial state (3c-I) and in the final state (3c-II) with respect to 
the hard interaction, in agreement with the general arguments in [E9j. 



11 



3 Gluon Contribution 



The gluon contribution to hard dijet production is more involved because of subtleties 
related to gauge invariance. One is tempted to calculate the coefficient function in front 
of the gluon distribution by considering pion scattering from on-shell gluons with zero 
transverse momentum, similar to the above calculation for the quarks. The result has in 
this case to be multiplied by the gluon distribution in the target in the physical light-cone 
gauge PiA^ = 0, which has the form 0, |23[ 



o 

where 



x x x 2 

(3.1) 



(pm) 

The difficulty arises because of the factor \/{x\X2) which is singular within the integration 
domain. The particular procedure to deal with this singularity is related to the gauge 
condition for the gluon field at time infinity and has to be established by considering 
carefully the light-cone limit of the relevant Feynman diagrams. For the classical process 
of hard vector meson production by longitudinally polarized photons || §] the correct 
prescription was formulated by Radyushkin 1 23[ 



X1X2 (x\ — ie)(x2 + ie) 



(3.3) 



It is not obvious, however, whether the same substitution is true for the hard dijet produc- 
tion. In fact, we will argue that no simple prescription exists in this case at all. In order 
to find the answer, in this section we will not assume zero gluon transverse momentum 
from the beginning, but consider instead the full scattering amplitude of the hard dijet 
production from the quark target mediated by the two-gluon exchange, cf. Fig. [[]. We will 
identify the (IR divergent) contribution to this amplitude corresponding to the region of 
small gluon transverse momenta and try to find a factorized expression for it, as a product 
of the coefficient function times the perturbative gluon distribution in a quark. 

For pedagogical reasons, we first consider the simpler case of hard vector meson pro- 
duction by longitudinally polarized photons. We do the calculation in Feynman and in the 
light-cone gauge and explain how the prescription in ( |3.3| ) arises in both cases. Next, we 
present the calculation of the imaginary part of the amplitude for hard dijet production in 
Feynman gauge. This calculation was previously reported by us in ||19|| . Finally, we derive 
the full expression for the amplitude (both real and imaginary part) using the axial gauge 
and compare our results to the work 



3.1 Getting Started: Exclusive Vector Meson Production 

To the lowest order in a s the relevant Feynman diagrams are shown in Fig. [7]. The 
kinematics is similar to that of the dijet production. We take P2 and p' 2 to be the target 
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Figure 7: The leading-order gluon contribution to the hard exclusive vector meson production 

* o 
7 ->P- 

quark momenta in the initial and the final state, respectively. The momentum of the 
vector meson in the final state is denoted by p\. The virtual photon momentum is equal 
to q = pi — (p 2 and we use the further notations Q 2 = —q 2 and ( = Q 2 /s with s = 1p\p 2 . 
We take the photon to be longitudinally polarized, with the polarization vector 



I 1 Q 1 Q L ii 

V = qPi,? + —P2,n = ^ + 2— p2,/i , 







4 



1 . 



(3.4) 



Thanks to the U(l) gauge invar iance one can omit the first term in the above expression 
and use a simpler vector 

(3.5) 



2—P2. 



The vector meson distribution amplitude is defined by the expression similar to ( |2.3| ) 



(0\q(0)^Q(y) \V(pi)) y *=o = fvfav) dze~ iz ^Mz) 



(3.6) 



where fy is the corresponding decay constant. For simplicity we do not elaborate on the 
isospin structure and consider below the contribution of the single light quark flavor q. The 
variable z corresponds to the momentum fraction carried by the quark. In the calculation 
of the amplitude sandwiched between the two target quark spinors u(p' 2 ) and u{j>2) we will 
assume averaging over the color of the target quark and make the replacement 



u(p 2 )M u(p 2 



Sp [M M 



(3.7) 



picking up contribution of the unpolarized gluon distribution, which is dominant at large 
energies. 

For definiteness, consider the contribution of the diagram in Fig. |a. We obt am 

l 

I f d^k f 

iM a = V4roe,3 4 /vVW^C f J -^—^ J dz<p v (z) 

o 



x Tr 



2Q 

s 



K Z Pl~ 



-ry 

[zpx - q) 2 + ie ((p 2 - k - zpi) 2 + ie 



7m2 



K 



(3i 
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where e q is the quark electric charge, a = 1/137 is the fine structure constant, and 



K 



2 



^27, 



N^ lUl (k)N^ U2 (k - (p 2 ) 
(pa - A;) 2 + it [k 2 + ie][(k - (p 2 ) 2 + ie] 



-7, 



(3.9) 



is a common factor for all Feynman diagrams in Fig. [j] that describes the emission of the 
pair of gluons from the target quark line. In the expression for TZ we have assumed a 
general form for the gluon propagator 



AT 

k 2 + ie V ' 



(3.10) 



In order to do the integration over the loop momentum we use the Sudakov parametriza- 
tion 

k = api + X1P2 + k± , d 4 k = -dadx\d 2 k± , k 2 = ax\s — k\ (3.11) 

and take the integral over a by residues. After this, the integration over x\ gets confined 
to the interval < x\ < 1 and the variable x\ acquires the meaning of the longitudinal 
momentum fraction. In the Feynman gauge N^ v = g^ v we obtain after some algebra 



M a = V4^e q f v ^/l - C 



Q C F g 4 



s 167nW r 



dz(J) V (z) / d 2 k± / dxi I a 



(3.12) 



with 



zs(2s(l — Xi) 2 z + k\x2) 



szx 2 - k 2 1 {l-C > )/{l-xx)+ie\ 

2s(l - xx)zx2 + k\(2 - xi - C) 



+<d(-x 2 



X2C 



(3.13) 



The contribution of large transverse momenta k± > Q in Eq. ( |3.12| ) gives rise to the 
a 2 correction to the coefficient function in front of the quark distribution in the target 
and is not relevant for our discussion. We concentrate, therefore, on the region of small 
k± <^ Q. Notice that the k±_ integral is infrared (IR) divergent. The quadratic divergence 
dk\/k\_ must cancel in the sum of all Feynman diagrams and the remaining logarithmically 
divergent integral ~ dk\jk\ has to be interpreted as the perturbative gluon distribution 
in a quark. In order to observe factorization one has, therefore, to expand the integrand 
in powers of k\ and keep two first terms. The expansion of the denominator in the first 
term in Eq. (|3.13| ) in the vicinity of x 2 = (x\ = () may seem dangerous. However, the 
integral over x\ can be deformed to the complex plane away from this singularity, and the 
expansion can be done without further problems. We obtain 

2s(l - x 1 ) 2 zx 2 + k\ ((1 - x x f + (1 - C) 2 ) 



+ e(-x 2 ) 



(x 2 + ie) 2 

2s(l - xx)zx 2 + k\(2 -xx-Q 
X2C 



O{k ± ) . (3.14) 
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Since the integration is restricted to the region k\ -C Q 2 the constant terms 0{kP L ) can 
be omitted. Note the +ie prescription to go around the singularity as a "memory" of 
the direction in which the analytical continuation to the complex x\ plane is performed, 
which, in turn, has its origin in the Feynman +ie prescription in the quark propagator 
1/[(CP2 -k- zpiY + ie}. 

Calculation of the other diagrams in Fig. [?| is similar. Here we cite only the final result, 

cf- 



M 



v 



where 



Am V '47T a e q a s f\ 

o 

l + (l-x) 2 -C 

i-C 



dz 



4>v{z) 



zz 



-e(c-x) 



(xi - ie)(x 2 + ie) 
(C-x)(2-x-Q 

C(i-C) 



dk] 



(3.15) 



(3.16) 



can be identified with the (perturbative) generalized gluon distribution of a quark|]. Owing 
to the +ie prescription to go around the singularity at x 2 = the amplitude acquires an 
imaginary part 



Im M 



-47r 2 V47ra e q a s fvs 



zz 



(3.17) 



Note that the l/kj_ terms and the double pole, ~ l/(x 2 + ie) 2 present in Eq. (|3.14j) cancel 
in the gauge invariant sum of Feynman diagrams, as expected. 

For comparison, let us do the same calculation in the light-cone gauge, with the prop- 
agator 



DUk) 



k 2 + 



(kpi 



(3.18) 



For the moment we do not specify a particular prescription to go around the singularity 
at kpi = 0, this choice will be discussed in detail in what follows. 

Using this expression we obtain the following result for the radiation factor (|3.9| ) in 
the light-cone gauge: 



K 



LC 



x x x 2 [k 2 + ie] [(k - (p 2 ) 2 + ie] [(p 2 - k) 2 + ie] ' 
where the numerator is equal to 



(3.19) 



7V%, = (2 - Qk^K + ]^asxix 2 g^ u + pi^p lu [Qak\/ s + 8a 2 (l - xx)\ 

+ k^p lv [2k 2 Js + a(4 - 3xi)] + k^p lt \2k\js + a(4 - 3xi - C)] (3.20) 

and the prefactor l/(xix 2 ) comes from l/(kp\) ■ l/[(k — CP2)Pi\ in the propagators 
of the t-channel gluons. 



§Eq. ( |3.15| ) is written assuming a single light quark flavor. For p° meson one has to substitute e q fv 
f p /y/2 with f p ~ 200 MeV. 
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The integration over the Sudakov variable a ( |3.11| ) converges at values a ~ k\/s. The 
numerator in ( |3.19| ) is, therefore, of order 0{k 2 1 ) and the radiation factor (|3.9| ) in the 
light-cone gauge has at most a 1/^ singularity at small transverse momenta, compared 
to l/kj_ in the Feynman gauge. Because of this, the QCD factorization in light-like gauge 
is valid for each diagram in Fig. ^ separately: In the upper parts of the diagrams one can 
substitute k = x±p 2 + api + k± — > x±p 2 and neglect the gluon virtuality and transverse 
momentum altogether. This is how the parton picture emerges: The amplitude for hard 
exclusive meson production is given by the convolution of the scattering amplitude off the 
on-shell transverse gluon 7*g — ■> Vg and the gluon distribution in the target quark. In 
perturbation theory the gluon distribution is given by the integral 



da 



d 2 k ± n]2- 



T 9 Ax x ) 



1 1 IV 



_c 

X1X2 



(3.21; 



In fact, only the two first terms in ( |jj.2(J ) contribute. Performing the integral over a and 
averaging over the directions of k± we indeed reproduce the result in Eq. 



The argument presented above is standard and tacitly assumes that the singularities at 
£2 — > and x\ — > play no role[]. In order to recover the correct prescription in Eq. (|3.15|) 
we have to be more careful. By the explicit calculation of the six diagrams in Fig. [7| we 
obtain the following expression: 



M 



where 



ia K 



d(as) dx\ d 2 k± [(2 — C)^l — (as)xiX2]y/l — ( 

k\ + ie] 



47T 3 J X1X2 [(as)xi — k\ + ie][(as)x2 — k\ + ie][(as)(zi — 1) 



• J (3.22) 



J 



27r\/47ra e q a s f\ 



dz (f>v{z) 



x 2 s 



x±s 



[(a + z)sx2 — kj_ + ie] [(a — z)sx\ — k\ + ie] 



[(a + z)sx 2 — kj_ + ie] [(a — z)sx\ — k\ + ie] 



+ (z *-> z) 



(3.23) 



The three terms in the square brackets in (|3.23| ) correspond to the diagrams in Fig. 0a,b,c, 
respectively, while the symmetric contributions (z <-> z) originate from the remaining three 
diagrams. If x±, \x2\ 3> k^/s then both a and k\ can be neglected in this expression and 
the factor J is easily recognized as the 7*^ — > Vg on-shell amplitude 



J 



dz 



<pv{z) 



zz 



(3.24) 



Note that in this case the diagrams in Fig. [7|c and Fig. |^f do not contribute. Performing 
the remaining integration and collecting all factors we recover the result in ( |3.15| ) obtained 
earlier in the Feynman gauge. On the other hand, the contribution of the singularity at 



^The singularity at x\ — > 
this point || , cf. Eq. ( ^16| ). 



is in fact irrelevant since the gluon distribution !Ff(x\) ([ 



y) vanishes at 
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x 2 = 0, or, equivalently, of the imaginary part of the amplitude requires some attention. 
First note that the expression for J vanishes identically when x\ — > or x 2 — > 0. As 
a consequence, the l/(x\X2) factor appearing in Eq. ( |3.22|) does not in fact produce any 
singularity in the gauge invariant sum of all Feynman diagrams. The result of the calcu- 
lation, therefore, does not depend on a particular prescription to go around the auxilary 
singularity in the gluon propagator in the light-cone gauge. Any prescription produces 
the same result. Since the poles at x 2 = and x\ — are spurious, we are left with the 
imaginary parts corresponding to +ie prescription in the Feynman propagator in (|3.23|) 
and invoking the standard argument with the contour deformation recover the x 2 + ie 
in ( p,15| ). The gauge prescription becomes important, however, if one insists on QCD 
factorization in each Feynman diagram separately. Consider x 2 — > and two different 
prescriptions for the factor l/[(k — CV2)Vi =t * e ] i n the gluon propagator. The first term 
in ( |3.23| ) corresponding to the diagram in Fig. 0a contains a factor x 2 in the numerator 



and is not affected. In the second term, corresponding to Fig. 0b, one can still take the 
light-cone limit and neglect a and k±. The contribution of this diagram then reads 

M [b) ~ / dx 1 , C , (3.25) 
J x 1 (x 2 ±ie)z 

It has both real and imaginary parts and the sign of the imaginary part depends on 
the gauge prescription. Note that in covariant gauges this 'crossed box' diagram has no 
imaginary part in the physical region. The third contribution in Eq. ( |3.23| ) originates 
from the diagram in Fig. 0c and is more delicate. This contribution vanishes in the 'naive' 
light-cone limit. In the vicinity of the point x 2 = the small factor k\ in the numerator 
is compensated, however, by the small denominator: 

Mm ~ dxx , c , '— = -. (3.26) 

J xi(x 2 ± ie) z [(a + z)x 2 s — k x + ie\ 

At small k\ and a ~ k\j ' s this reduces to 

T\{x x ) 1 ( 1 



M „ / dx 1 ^—- - r — T + r — tt • (3.27) 
J xi z V [x 2 ± ie\ [x 2 + ie\ J 

We see that first, the gauge poles [x 2 ±ie] indeed cancel in the sum of diagrams in Fig. 0b 
and Fig. 0c and, second, the diagram in Fig. 0c is equal to zero in the light-cone limit 
if and only if one chooses the +ie prescription for the gauge pole. The reason becomes 
obvious if one tries to deform the integration contour in x± away from the singularity at 
x\ = ((x 2 = 0), as shown in Fig. ||. If the +ie prescription for the gauge pole is used, both 
the gauge pole and the pole of the quark propagator lie below the real axis. Hence the 
contour deformation is not obstructed and taking the limit k± —>■ we obtain zero, the 
parton model result for the contribution of this diagram. On the other hand, if the — ie 
prescription is used, the gauge pole appears to be above the real axis and in the k± — » 
limit is pinched with the pole of the quark propagator. In this case one cannot move 
the contour away from the singularity. As the result, the diagram in Fig. 0c acquires a 
non-zero imaginary part which has no parton model interpretation and, most importantly, 
is missed in the 'naive' calculation when k± and a are put to zero at the beginning. 
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The net outcome of our discussion is the verification of the ansatz in Eq. ( |3.3| ) for the 
leading-order contribution in the strong coupling: One can calculate the hard coefficient 
function as on-shell 7' *g — > Vg amplitude, and use the prescription in ( p.3|) for the definition 
of the gluon distribution in Eq. (|3.1|) . In order to assemble the all-order proof, one has 
to show that to arbitrary order in perturbation theory the poles of Feynman propagators 
corresponding to soft gluon exchanges all lie below the x 2 real axis. The corresponding 
discussion can be found in |7[]. The physical interpretation is that all soft gluon exchanges 
can be thought of as part of the final state interaction between the outgoing hadrons; 
they are reduced to eikonal factors that enter the (gauge- invariant) definition of the gluon 
distribution. Pinching of the singularities in this case cannot occur since initial state soft 
interactions are not present. The interpretation of the ±ie prescriptions in ( |3.3| ) as due to 
soft final state interactions allows for an alternative derivation of this result, by noticing 
that the suitable boundary condition for the gluon field in the light-cone gauge is such 



that [31, 23 



Mv) =PiJ G AV + <ypi)e- ea da (3.28) 


with the integration extended to plus infinity in time. We remind that pi is the pion 



momentum. It is easy to check |23] that the ansatz in ( |3.3| ) is an immediate consequence 
of this relation. A still another interpretation of this ansatz is that the amplitude of hard 
exclusive meson production only has a s-channel discontinuity in the physical region. Its 
energy dependence has, therefore, to be a function of s + ie which translates to £ — > £ — ie 
and x 2 — > x\ — ( + ie = x 2 + ie, respectively. In any case, we see that the ansatz in ( ^.3| ) is 
specific for the considered process and cannot be taken over for the hard dijet production 
without a careful analysis. In fact, the very existence of such an ansatz is a consequence 
of QCD factorization^} 

3.2 Dijet production: Dispersion approach 

Since the theoretical status of hard dijet production continues to be controversial |T7|, 



18] , [19], PD| , in this paper we will present our calculation using two different techniques 
and show that both lead to the same result. In this section we calculate the imaginary 
part of the 7r — > 2 jets scattering amplitude on a quark target in the high energy limit. 
The imaginary part is interesting in several respects. First, we have found in Sec. 2 
that the imaginary part of the quark exchange is affected by the pinch singularities and 
contains logarithmic end-point divergencies which destroy collinear factorization. Since 



II An instructive example is provided by the process Vp — > "f*{Q 2 )p where the photon in the final state 
has positive virtuality. This reaction is similar to electroproduction and can be treated in the same way. 
In this case, however, the leading-order calculation yields the prescription [xa — ie], i.e. opposite to the 



ansatz in (3.3). The reason is that the soft interaction is in the initial state in this case, and in order 
to preserve the parton picture one has to use a different gauge condition. Another difference is that the 
amplitude Vp — > 7*(Q 2 )p has nonzero dispersive parts in both variables s and Q 2 . One can check that 
the contributions of s-channcl singularities all cancel in the light-cone limit with part of the singularities 
in the Q 2 -channel. Hence one is left with Q 2 -channel singularities alone, and the dependence of the 
amplitude on [x\ — ( — ie] can be understood as the replacement Q 2 — > Q 2 + ie. 
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Feynman diagrams describing quark and gluon contributions have similar topologies, the 
same problem is expected for the gluon contribution as well. Second, at high energies the 
scattering amplitude corresponding to Pomeron exchange is dominated by its imaginary 
part, so that this contribution is numerically the most important one. Last but not least, 
the cut diagrams appearing in the calculation of the imaginary part are built of tree-level 
on-shell scattering amplitudes and their form is strongly constrained by gauge invariance, 
see below. This simplification has been widely used in the literature in the calculations of 
high-energy asymptotics of scattering amplitudes starting from |32|, [33j. The idea to use 
this approach for the dijet production was suggested in |L8 . 

The s-channel discontinuty of the amplitude irq — > 2 jets + q is described by the cut 
diagrams shown in Fig. ^. They can be grouped into the four gauge-invariant contribu- 
tions in Fig. Rla-d which differ by the position of the hard gluon that provides the large 
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Figure 8: The decomposition of the imaginary part of the amplitude nq 
gauge-invariant contributions. 



(qq)q into four 



momentum transfer to the jets. The corresponding contributions to the imaginary part 
of the amplitude will be denoted by Cr a ), Co,), C( c \ and C^), respectively. For example, in 
Fig. |a it is assumed that the hard gluon exchange appears to the left of the cut. This 
contribution is given by the sum of 10 Feynman diagrams, for further details see our letter 
|T9|| . Similarly, the contribution in Fig. |8|b is given by the sum of 10 diagrams with the 
hard gluon exchange appearing to the right of the cut. The two remaining contributions 
in Fig. |8]c and Fig. |8]d take into account the possibility of real gluon emission in the inter- 
mediate state. The filled circles stand for the effective vertices describing gluon radiation, 
see Fig. |9|. Each of the two contributions in Fig. ||c,d corresponds to a sum of 9 different 
Feynman diagrams. 




Figure 9: The effective vertex. 



At first sight, we have to consider the discontinuity in the invariant mass M 2 of the 
jets as well. It is easy to check that the corresponding cut diagrams have to have the hard 
gluon exchange to the right of the cut in order that the transition from the intermediate 
to the final state corresponds to a physical process. On the other hand, unlike for the 
s- discontinuity, the t-channel gluon lines can be crossed in this case. The two possible cut 
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b") b") 
Figure 10: The M 2 -discontinuity of the amplitude irq — > (qq)q. 

diagrams are shown in Fig. [10]. We denote the corresponding contributions by Ciy) and 
C(6«), respectively. They differ by the interchange of the Mandelstam variables s and u. 

It is known for a long time that at high-energies any amplitude with vacuum quantum 
numbers exchange in the t-channel has to be crossing-symmetric, or in other words has 
positive signature. In our case we are dealing with a two-gluon exchange that grows 
linearly with energy and consider a one-loop amplitude which adds a logarithm. The 
crossing symmetry then implies that the contribution to the scattering amplitude of which 
C(b), C(w) and Cn,") are the three nonvanishing dispersion parts in the physical region of the 
s-channel has the following schematic structure: 

M ~ s In + u log (-^) + O (s°) . (3.29) 

At high energies u ~ — s and the above structure implies that that the s-channel and 
the M 2 -channel discontinuities are related: C^) = — C(b') and C^>) = —£(&")■ The first 
equality = — C(ft') which amounts to the cancellation of the s- channel and M 2 - channel 
discontinuities is not specific for high energies, but rather a general property of the scaling 
limit: In so far as the amplitude is a function of the dimensionless ratio (— s)/(— M 2 ), 
the s- and M 2 - discontinuities cancel each other identically^. It follows that the total 
imaginary part of the amplitude 7tq — > (qq)q in the scaling limit is given by the sum of the 
three s-channel cut diagrams in Fig. |]a,c,d and the M 2 -cut diagram with crossed t-channel 
gluon lines in Fig. |T0|b". However, at high energies the crossing symmetry implies that 
C(y/) = C(b) and the net result is that the M 2 -discontinuity can be neglected altogether. 
We are left, therefore, with the set of cut diagrams shown in Fig. ^. 

The general strategy of the calculation is the following. As is easily checked by inspec- 
tion, in any cut diagram two internal lines are on the mass shell. The corresponding two 
on-shellness conditions fix a and x\ variables in the Sudakov parametrisation ( |3.11| ) of the 
gluon momentum and relate the variables z' and X\ (see Fig. [I] for the notations) to one 
other. Since the a and x- variables for both gluons are of the order of 1/s, the 1/& 2 2 factors 
in the propagators of the t— channel gluons can be approximated by using k 2 ~ — fc 2 , to 
the 0(1/ s) accuracy. In this calculation we use the Feynman gauge and perform the usual 
substitution g^ v — > P2V1HP1P2) in the numerators of the t— channel gluon propagators, 
which is valid to the same accuracy. Using the on-shellness conditions for the contribu- 
tions in Fig. |8]a and Fig. ^0 one obtains x\ = (, X2 = 0, for any z' . For Fig. |S]d one finds 
x\ = (z'z/(z' — z), X2 = (zz'/(z' — z) and z' > z, where the last condition ensures that 

**We have already met with an example of such cancellation for the quark contribution, see the discus- 
sion of cut diagrams in Fig. ^(2c-I) and Fig. ^(2c-II). 
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the energy of the cut gluon is positive. Finally, for the set of cut- diagrams corresponding 
to Fig. |8|c we obtain x\ = (zz 1 / (z — z'), x 2 = (z'z/ (z — z') and z > z' . 

After the on-shellness conditions are used, a single integration over the gluons trans- 
verse momentum k± remains: 

/d 2 k± 
jpy Ju P (k±, q±) Jdown(k_L, q±) , (3.30) 

where k\ comes from the product of the two gluon propagators. J up and J down are dubbed 
impact factors and stand for the upper and the lower parts of the diagrams in Fig. |8|a-d 
that are connected by the two-gluon exchange. The representation Q3.30| ) is well known 
|32| , |33| from studies of QED scattering at high energies. 

Properties of the impact-factors J up and Jdown as a functions of k±_ at k± — > are 
of crucial importance. Since Jdown is the impact-factor of a point-like target quark, 
Jdown(kj_, q±) ~ const. On the other hand, J up {k±,q±) stands for the scattering of the 
colorless qq (Fig. |]a-b) or qqG (Fig. ||c-d) state having a transverse size ~ l/q± and has 
to vanish at small fci < gi, J U p{k±, q±) ~ kj_, as a consequence of the color neutrality of 
the quark-antiquark pair: A gluon with a large wave length ~ l/k± cannot resolve a color 
dipole of the small size ~ l/q±- Since in our case there are two gluons, J up is proportional 
to the product k±-k± = k^j. In the opposite limit of large transferred momenta, k± ^> q±, 
the two t— channel gluons are forced to couple to the same parton (quark or gluon) in the 
upper block in Fig. [8]a-d. It follows that at large k± J up (k±, q±) ~ const. 

Taking into account the properties of the impact-factors discussed above, we conclude 



that the transverse momentum integration in (|3.30| ) diverges logarithmically at small k± 

2 

and the integral can be estimated by Ai ~ f q± dk\jk\ ~ lng^, as expected. The region of 
k\ > q\ does not produce the large logarithm and can be neglected. Note that the correct 
small k± behavior of the impact factors is only recovered in the sum of cut diagrams for 
the gauge invariant amplitudes C( ), Cq,), C(c) and C^, but not for each diagram separately. 

In addition to the diagrams discussed so far, the amplitude nq — > (qq)q receives a 
contribution from the three-gluon exchange in the t-channel. Such terms can be viewed 
as belonging to the cut diagrams shown in Fig. [8]a in which the hard gluon in the blob is 
attached to the bottom quark line. We have checked that this extra contribution does not 
contain the large collinear logarithm lnq\ and therefore we neglect it. 

The calculation of C^) proceeds as follows. Let P = aip 1 ^ + xip% + 1± be the momentum 
of the (real) gluon in the intermediate state and let e^(l) be one of the two physical 
polarization vectors. The two conditions (e-p2) = and (e- 1) = fix the gauge and result 
in e"(Z) = e£ + 2p£ (ej ± ) /(a, s). 

The effective vertex corresponding to the sum of the three diagrams in Fig. |9| has the 
form 



.g 2 z(z' -z) 



I ( t H a )..-- (t a t 1 ).. 



u(qi 



z' — z 



V -u{z'p x ). (3.31) 

s 



q z ± z' 

Here t l and t a are the SU(3) generators. The color indices I and a belong to the emitted 
^ The 0{k\) behavior can be traced to the gauge invariance of the amplitude, see p3] for the details. 
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gluon and the t— channel gluon, respectively. We have also introduced an auxiliary two- 
dimensional vector defined as: 



6* = A £_ 2 (^) tf 
q{ 



ft. 



(3.32) 



Note that the effective vertex, in the limit of small k±, is proportional to b oc k±. The 
constant terms cancel in the gauge invariant sum of the diagrams in Fig. ||. 

The second effective vertex in Fig. |8]d has a similar form. Combining both of them 
and performing the sum over the polarizations of the emitted gluon we obtain the impact- 
factor jip . Since each effective vertex is proportional to k±, it follows that J$ ~ kj_, as 
expected. The result for the amplitude C(d) is obtained using the representation in fl3.30p . 
The calculation of C( c ) is very similar. The result for their sum reads: 



C( c ) + C 
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zz 



DC 2 F 



1 + 



dk\ 
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dz' 4> n (z') 



(N- 



z z 



+ 1 x 



Q{z' - z) Q(z - z') 
{P - z) + (z - V) 



(3.33) 



where 



D 



■ f 3 4tt 2 
~ lSUas N?qi 



u{qi) r r 5 —v{q2) 5ij8 c 



(3.34) 



The color indices correspond to the produced quark-antiquark pair (jets) and (c, c') 
stand for the color indices of the target quark in the initial and the final state. The 
contributions ~ Q(z' — z) and ~ Q(z — z') belong to M.(d) and -M( c ), respectively. 
For the cut diagrams in Fig. Ha and Fig. Rib we obtain: 
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DCl 



dk\ 
*i 

dk\ 
~~kT 
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dz'(f> n (z') [j f + j t) . 



dz 



z z\ 1 f zz 

:: | i 7 + i 7 J + (iV 2 -i) l^i 7 



+ 1 



. (3.35) 



The transverse momentum integrals in ( |3.33| ) and ( |3.35| ) can be identified with the 
small x limit, i,( < 1, of the (perturbative) generalized gluon distributions of a quark 
defined in Eq. ( |3.16| ): J-f!(x) ~ ^Cp f dk\jk\. We obtain, in this approximation |I9[ 



-Im [iM g i uon ] 

7T 



s f* a l 4^~T M (?i)75— v(q 2 ) TS, 
JS„ q i s 
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(3.36) 



with 
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0(z' - z) 



(z'z 



0(z 



T9 



(z'z 



(3.37) 



The expression in Eq. Q3.37| ) presents the main result of this section. Note that the 
integrand in ( 3.37 ) is singular at z' = z so that there is a logarithmic enhancement of 
the contribution of the integration region £ <C \z' — z\ <C 1. In addition, there is a 
logarithmic divergence at the end-points z' — > and z' — > 1 which signals that the collinear 
factorization is broken, as expected. In what follows we will discuss the contributions from 
these regions in some detail. Before doing this, however, in the next section we derive the 
complete result for the amplitude in the scaling limit (both real and imaginary parts, and 
including C(l/s) corrections) using a different approach. 



3.3 Dijet production: Factorization and the light-cone limit 

Our aim in this section is to derive the complete result for the leading-order contribution 
to the kernel Tfj(z, z r , x±, x 2 ) such that 

1 1 

A-K 2 a 2 s sif n n — N _ i> 2 



-M g i U on = — N * 4 * V 1 ~C^(gi)75y^(g2)^y dz' (p^z') J dx 1 J 7 %(x 1 )T gluon (z,z',x 1 ,x 2 ), 
c ± 00 

(3.38) 

cf. Eq. (|1 . 2| ) . For high energies, £ — > 0, we expect to recover in this way the result of the 
direct calculation of the imaginary part given in ( |3.37| ). 

We will mainly be concerned with the singularity structure of T g i uon (z, z' , x%, x 2 ) at 
X2 — > 0. The singularity structure at x\ — > can, in principle, be found from similar 
considerations. This is in fact not necessary since it can be established using the crossing 
symmetry T g \ uon (z, z', x 1; x 2 ) = T g \ uon (z, z', —x 2 , —xi) that corresponds to the interchange 
of s- and w-channels for the corresponding pion-gluon amplitude. Apart from this issue, 
the calculation of T gluon (z, z', xi, x 2 ) is straightforward and can most easily be done by 
considering pion scattering from on-shell transversely polarized gluons, cf. EQ]. For the 
simpler case of hard exclusive production of vector mesons we have argued in Sect. 3.1 that 
the light-cone calculation (with on-shell gluons) is actually sufficient since the singularity 
structure in Eq. ([3.15 ) can be restored using the prescription ( |3.3|) in the definition of the 



off-forward gluon distribution (|3.1p. One way to understand this result was that for hard 



exclusive production of vector mesons it is possible to choose an axial (light-cone) gauge 
in such a way that gauge singularities of the t-channel gluons and causal singularities of 
Feynman propagators lie on the same side of the integration contour in x% so that there 
are no pinches in the Glauber region. In the present case, a similar simplification is not 
expected since soft gluon exchanges occur both in the initial and in the final state. 

The complete set of relevant Feynman diagrams is shown in Fig. [TT|. For 20 out of in 
total 31 existing diagrams, Fig. |TI](l-20), it is kinematically possible that the t-channel 
gluon with vanishing momentum fraction x 2 couples either to the initial or to the final 
on-shell quark lines. The corresponding quark propagators produce a singularity in the 
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x 2 — > limit. In particular, in the 5 diagrams in Fig. 0(11-15) where the soft gluon is 
attached to the outgoing quark, the quark propagator is 

l/[(?i + k-(p 2 ) 2 + ie] = l/[(z + a)(x 1 -(z)s-(q ± + k ± ) 2 +ie] -> l/[zs(x 2 + ie)} (3.39) 

in the limit a, k± — > 0. Similarly, in the 5 diagrams in Fig. |Tl|( 16-20) where the gluon is 
attached to the outgoing antiquark, one obtains 

l/[{q2 + k-(p 2 ) 2 + = l/[(z + a)(x 1 -Cz)s-(q ± + k ± ) 2 + ie] -> l/[zs(x 2 + ie)} (3.40) 

in the same limit. However, in the 5 diagrams in Fig. |Il](l-5) where the soft gluon is 
emitted by the ingoing antiquark, the denominator of the quark propagator is 

l/[(z' Pl — k + (p 2 ) 2 + ie] = l/[(z' - a)(C - x 1 )s -k\ + ie] -> l/[-z's(x 2 - ie)} (3.41) 

and similar for the diagrams in in Fig. |Tl"|(6-10) 

l/[(z' Pl - k + (p 2 ) 2 + ie] = l/[(z' - a)(C - x ± )s -k 2 ± + ie] -> l/[-z's(x 2 - ie)] . (3.42) 

We see that both types of singularities are produced, (x 2 ± ie) , depending on whether the 
soft interaction occurs in the initial or in the final state [|29|| . As a consequence, it is not 
possible to choose an axial gauge in which the light-cone limit exists for each Feynman 
diagram separately. 

For definiteness, let us define the light-cone gauge in the same way as this was conve- 
nient for hard exclusive production of vector mesons, l/x 2 — > 1/ (x 2 + ie) . In this case the 
diagrams with soft interaction in the final state Fig. [TT|(1 1-20) and the diagrams that are 
non-singular in the x 2 — > limit Fig. 0(21-31) can be calculated in the light-cone limit, 
with on-shell t-channel gluons, and the result multiplied by the gluon distribution ( |3.1| ) 
using the prescription in ( p.3| ). On the other hand, in the diagrams with soft interaction 
in the initial state Fig. |Tl|(l-10) the integration contour over x± will be pinched. We can 
calculate their contribution by deforming the integration contour in the same direction as 
for the other diagrams x 2 = x\ — ( + ie, and adding the contribution of the pole of the 
quark propagator fl3.41| ) and ( |3.42| ). In what follows we refer to these two contributions 
as 'naive' light-cone T^ LC and pole T? ole contributions, respectively. The subscript 'i' nu- 
merates the diagrams according to Fig. [11]. The additional pole contributions are present 
(with our gauge choice) in the diagrams in Fig. |TT](1-10) so that 

31 10 

T gluon (^, z', x u x 2 ) = T LC + T? ole = T t° + E T ?° Xe ■ ( 3 ' 43 ) 

i=l i=l 

T LC can be calculated as the gauge-invariant on-shell pion-gluon scattering amplitude, 
times the gauge factor ( p.3|) 

t lc = £ y x ^m 9 (3>44) 

{xx -ie)(x 2 + ie) ^ 



where A4 correspond to the upper parts of the diagrams in Fig. TT], not including the target 
quark and the propagators of the t-chanel gluons. Since ^9-^{m)g = J^ 3 ^ l 4j r9 " > ™ is the 
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amplitude of a physical process, it is gauge invariant and can be calculated in any gauge. 
In the Feynman gauge we find 
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The remaining contributions are given by the replacements 



(3.45) 
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(3.46) 



For the sum of all diagrams we obtain 



T LC (z,z',Xx,x 2 ) = C F (^ + j) 1 C 

.z'z' J 
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[xx - ie] 2 [x 2 + ie] 2 [x\ - ie] [x 2 + ie] 
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(3.47) 



ZZ ( z z 



2N c z'z' \z'z' 



zz 



[xx - ie] [x 2 + ie] 



Taking the imaginary part of Q3.47| ) we reproduce the answer for the imaginary part of 
the coefficient function obtained in Ref. |20Jq 

In order to make the comparison it is necessary to take into account that the answer given in |f2of is 
rewritten in a different form using the symmetry (f> T (z) = 4> v (l — z) of the pion distribution amplitude. 
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As explained above, this is not the whole story, however, and one has to add the pole 
contributions. For example, consider the diagram in Fig. |lT](3) . In the notation of Sect. 3.1 
we obtain the corresponding contribution as 



M, 



c\ 



d 4 k 



dz' Mz')KL N ^( k + z Vl - q x ) (3.48) 



(2tt)* 
o 



(p 2 - z'pi) 2 + ie 



where 7^uf„ 2 is defined in Eq. (|3.19| ) and N ulV2 {k + z'pi — qx) is the nominator of the hard 
gluon propagator in axial gauge. In the present case only the first two terms in Eq. (|3.20|) 
contribute and according to our procedure we choose the prescription l/x 2 — > l/[x 2 + ie] 
in the corresponding factor in Eq. (|3.19|) . This singularity is due to gauge fixing and is 
pinched with the quark propagator [(k — (p 2 — z'pi) 2 + cf. ( |3.41| ). Although at the end 
we have to isolate the quark propagator pole contribution, it is convenient to rewrite the 
gauge pole as 

' ' (2iri)6(x 2 ) (3.49) 



[x 2 + ie\ [x 2 

and find the corresponding two contributions to ( |3.48|) separately 



M 3 = M% + M 



(3.50) 



as an intermediate step. 

Calculation of the first term, .M3, is straightforward since poles of the quark propagator 
and l/[x 2 —ie] do not pinch. One can neglect k± and a in the integrand of ( 3.48j ) everywhere 
except in 7£^p . The integral of 7£„fL over a and k± reproduces the generalized gluon 
distribution, cf. ( |3.21| ), and gets factored out. The remainder gives the contribution to 
the coefficient function (in the light-cone gauge) 



T a = Cr Z^eig2z> - Z) - Xl ] 

z'z'xi[x 2 — ie\[(z' — z)x\ — (z'z + ie\(x\ — (z) 



(3.51) 



The factor (x\ — (z) in the denominator of the above expression originates from the hard 
gluon exchange propagator in the light cone gauge, the factor N UlU2 (k + z'pi — q±) in 
Eq. ( gggp . 



In the second term, A4\, we use the 5-function to put x\ — > ( and need to expand 



the quark propagators in Eq. (|3.48|) in the limit k± — > and a ~ k\/s. To the required 
accuracy 
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(3.52) 
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The integration over a produces the expression which can be identified with the pertur- 
bative leading-order generalized gluon distribution of a quark at the point J 79 {x\ = () , cf . 
Eq. ( |3.16|) . The corresponding coefficient function takes the form 



It 



(2m)S(x 2 )C 1 



z 3 



(3.53) 



As the final step, we go from l/[x2 — ie] in eq. Q3.51|) to l/[x2 + ie] using ( p.49| ) . The result 
equals the contribution of this diagram in the 'naive' light-cone limit, T% c , and the extra 
terms ~ S(x2) correspond to the 'pole' contribution T^ ole in the sense of Eq. ( |3.43| ). We 
obtain 
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(3.55) 



Note that in difference to the expressions collected in (|3.45|) the result in (|3.54j ) is obtained 
in the light-cone gauge. 

Calculation of the other diagrams is similar. We obtain 
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(3.56) 



and 
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(3.57) 



In these expressions 
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(3.58) 



The corresponding contributions originate from the last term in Eq. ( |3.20| ) and at first 
sight are dangerous because an additional factor a in the numerator of (|3.58|) makes the 
integral divergent at a — > oo. This difficulty is, however, spurious since the terms ~ rj 
actually cancel in the sum of all ten diagrams. For the sum of pole terms we obtain 
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and the final answer for the gluon coefficient function ( |3. 3S| ) can be written as 
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(3.60) 
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We see that, as we already discussed before, in this expression there are terms with different 
prescriptions for the poles at X\ = 0, x 2 = 0. In particular, in comparison to the light-cone 
result in ( |3.47|) , there are different ±ie prescriptions in the last line. 
The corresponding imaginary part is 
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(3.61) 



In the approximation that only the (large) imaginary part of the amplitude is taken into 
account, the differential cross section summed over the polarizations and the color of quark 
jets is given by 



dq 2 ± dtdz 8N^q 8 ± 



(3.62) 



The factorization scale /i 2 has to be of order of the transverse momentum of the exchanged 
gluon. 

The expressions in ( 3.60Q and ( |3.61|) present the main result of this paper. The ex- 
pression for the imaginary part in (|3.61|) indeed reproduces the result in (|3.37 ) apart from 
the extra last term containing a derivative of the gluon distribution. Using the explicit 
expression for the perturbative gluon distribution function of a quark in Eq. ( |3.16| ) it is 
easy to check that the derivative term is of order 0(C) and this is the reason why it was 
absent in the first calculation. 
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Note also that T pole given in Eq. ( |3.59| ) is nothing else but the M 2 -channel discontinuity 
C(hn\ given by the crossed diagram in Fig. [10 : 
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DC F — dz'<j) n {z') / dxiF?(xi 



-ImT polc 
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(3.63) 



cf. Eq. ( |3.35| ), up to 0{Q corrections. The imaginary part corresponding to the 'naive' 
light-cone limit is given therefore by the expression ~ ^(<j) — ^(&) ^( c ) 
reversing the sign of the contribution in Fig. (||)b. 



3.4 Discussion 



The leading end-point behavior of the coefficient function T g i uon (z, z', x±, X2) ( |3.60| ) at 
z' — > and z' — ► 1 is given by the following expressions: 

Tgl uon |^o -> "(2«)^ ^) [^l) +^2)] , 

T g i uon | z ^i -> -(27ri)|| (c F z + ^) [*(an) + 5{x 2 )\ . (3.64) 

Since T®{x) vanishes at x = [p3|l , only the terms ~ ^(^2) contribute. Note that the 
leading 1/z' 2 , 1/z' 2 asymptotics is imaginary, it comes entirely from the term which we 
denoted T pole , see Eq. ( |3.59| ). The contribution of T LC is at most 0(l/z'), or 0(l/z'), 



see also pOfl . Eqs. ( |3.64| ) have to be compared with Eqs. ( |2.13| ),( |2TT4] ) for the quark 
contribution. Similar to the latter case, the collinear factorization is violated due to the 
pinching of M 2 -channel and s-channel singularities in the Glauber region. Summing the 
both contributions in Eq. ( |3.64p and using the symmetry of the pion distribution amplitude 
we obtain [1115 1 



1 



end— points \ -''c/ J Z 



N c + ^)zz I dz' ™±jJL±J*(C, ^) . (3.65) 



Since (j)^{z') ~ z' at z' — > 0, the integral over z' diverges logarithmically. Remarkably, the 
integral containing the pion distribution amplitude does not involve any ^-dependence. 
Therefore, the longitudinal momentum distribution of the jets in the nonfactorizable con- 
tribution is calculable and, as it turns out, has the shape of the asymptotic pion distri- 
bution amplitude <fif{z) = 6zz. The corresponding physical process is the following. The 
limit z' — > 1 corresponds to a kinematics in which one quark carries the entire momentum 
of the pion. The fast quark radiates a hard gluon which carries the fraction (1— z) of quark 
momentum. This radiation is perturbative and is described by the effective vertex (|3.31|) 
at z' — 1. In the final step the hard gluon transfers its entire longitudinal and transverse 
momentum to the quark jet, and emits a soft antiquark which interacts nonperturbatively 
with the target proton and the pion remnant. It follows that the divergent logarithm in 
J 07r(-2') / z' 2 is of the form In q\/ /i 2 R where /ipR is related to the average transverse momen- 
tum of the quarks inside the pion. It is possible that in the case of scattering from a heavy 
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nucleus /zir may grow as ~ because of color filtering. On the other hand, it is likely 



that the modified factorization approach of [^6] can be developed and applied to coherent 



diffraction as well. A detailed discussion of these possibilities goes beyond the tasks of 
this paper. The other important integration region is the one when ( <C \z' — z\ <C 1, i.e. 
when the longitudinal momentum fraction carried by the quark is close (for high energies) 
to that of the quark jet in the final state. The enhancement of this region comes from the 
(small) denominator l/(z' — z) which is present in the contributions in Fig. ^p,d with real 
gluon emission in the intermediate state. The logarithmic integral J dz'/\z — z'\ ~ Ins is 
nothing but the usual energy logarithm that accompanies each extra gluon in the gluon 
ladder. Its appearance is due to the fact the the gluon in Fig. |]c,d can be emitted in a 
broad rapidity interval and is not constrained to the pion fragmentation region. To loga- 
rithmic accuracy we can simplify the integrand in ( |3.37| ) by assuming z' — z everywhere 
except for the small denominators and the argument of the gluon distribution, to get 



C<|z'-z|<l 



1 1 
AN e Mz) f-^-^-T-^ ql)^AN c Uz)l-H^<ll)- ( 3 ' 66 ) 

j z' — z ^ z' — z J y 



For a flat gluon distribution T^{y) ~ const at y — > 0, and the integration gives const -In l/( 
which is the above mentioned energy logarithm. Note that the color factors combine to 
produce Ca = N c signaling that the relevant Feynman diagrams in Fig. |8|c,d are those with 
a three-gluon coupling. Moreover, the factor 2N c /y appearing in ( |3.66|) can be interpreted 
as the relevant limit of the DGLAP gluon splitting function 



i i 



ql W^ {X = C " ql) = ^ I ^ P f(C^)^'^)^^ / dy^^(y,ql). (3.67) 



c c 



The quantity on the l.h.s. of ( |3.67| ) defines what can be called the unintegrated generalized 
gluon distribution and the physical meaning of Eqs. (|3.66|) and (|3.67| ) is that in the region 



z' ~ z hard gluon exchange can be viewed as a large transverse momentum part of the 



gluon distribution in the proton, cf. [17]. This contribution is proportional to the pion 



distribution amplitude 0„-(z,g^) and contains the enhancement factor ~ In s/qj_. 



4 Numerical Analysis 

Since collinear factorization does not hold, the numerical analysis presented below has to 
be considered as semiquantitative. We will use the simplest prescription to regularize the 
infrared divergence of the end-point contribution by imposing an explicit cutoff on the 
quark momentum fraction in the pion 

i i-MmM 
dz J dz (4.1) 
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Figure 12: Different models for the pion distribution amplitude, see Eq. Q4.2|), 



with /zir = 500 MeV. We use parametrizations of the generalized quark and gluon parton 
distributions by Freund and McDermott [B5| that are based on the new MRST2001 leading- 



order forward distributions |36j , see Appendix A, and the following trial pion distribution 
amplitudes: 

lCZ 



2 ( z ) = <fc L (z^ = 0.5 GeV) = 302(1 - z)(l - 2z) d 



3 ( z ) = (j>^{z, fx = 2 GeV) = 15z(l - z) [0.20 + (1 - 2z) 



(4.2) 



see Fig. [12]. The first expression in (|4.2|) is the asymptotic distribution amplitude at 
large scales H, the second is the Chernyak-Zhitnitsky model [3~4"| and the third is the 
Chernyak-Zhitnitsky model evolved to the typical scales fi ~ q± = 2 GeV probed by the 
E791 experiment f!4], [I5|. (j)i(z) and 4>?,{z) have very different shape and can be considered 
as two extreme cases. Our aim is to find out whether this difference can give rise to sizeable 
effects in the momentum fraction distribution of the jets. 

Different contributions to imaginary and real parts of the amplitude of the hard dijet 
production at s — 1000 GeV 2 and q± = 2 GeV from the nucleon target for the two 
extreme choices of the pion distribution amplitudes are plotted in Fig. [13] and Fig. \A . 
The normalization corresponds to Eqs. ( P-5|) , ( P-6|) , see Appendix B where we collect all 
the necessary expressions. It is seen that both the imaginary and the real part of the 
amplitude are dominated by the contribution of the generalized gluon distribution. The 
singlet quark contribution is less than 25%, and the contributions of w-quark annihilation 
and the <i-quark exchange are negligible. In |2TJ the quark contribution was found to 
be as important as the gluon one. One reason for the difference is that we use a more 
realistic parametrization for the quark distributions in our analysis. The second, more 
important reason, is that our result for the imaginary part of the gluon contribution 
contains an additional non-factorizable term, the last line in Eq. (|B.8| ), that is related 
with the discontinuity of the amplitude in the M 2 -channel. This additional contribution 
is large and constitutes approximately one half of the total result for ImJ gluon . Note also 
that the momentum fraction dependence of Im J g i uon deviates considerably from the shape 
of the pion distribution amplitude. This means that the large-energy approximation in 
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Figure 13: Different contributions to imaginary (left) and real (right) parts of the amplitude 
of the hard dijet production at s = 1000 GeV 2 and q± = 2 GeV from the proton target, as- 
suming the asymptotic pion distribution amplitude ^(z) = <fti(z), Eq. ( f4.2[ ). The normalization 
corresponds to Eqs. ( |B.5| ), (|B.6| ). Shown are: the gluon contribution (thin solid curve), the quark- 
singlet contribution (dashed curve) and the u and d-quark contributions (dash-dotted curve) . In 
addition, the sum of all contributions is shown by the thick solid curve. 
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Eq. ( |3.66| ), where <j> w (z) enters as an overall factor, is not adequate in the energy range of 



the E791 experiment. 

The real part of the amplitude is significant. The ratio Re/Im is not much smaller 
than unity and depends both on the assumed shape of the pion distribution amplitude 
and the jet momentum fraction. Qualitatively, the Re/Im ratio can be estimated using 
the Gribov-Migdal formula |37] 



ReM(s) 7i d , ImM(s) .. oX 

— = In — . (4.3) 

ImM(s) 2 dins s 

In our case, the amplitude is roughly proportional to the value of the generalized gluon 
distribution at the point X\ = (, M. ~ C _1 ^(C)> an d in the relevant kinematic range 

^(C)^0.52r - 43 , (4.4) 

cf. Appendix A. This gives ReM/ImAi ~ 0.43 • tt/2 ~ 2/3, in rough agreement with 
Figs. |13|Jl"4} The large real part is a yet another indication that the energy of E791 
experiment is not sufficiently high to consider the process as a diffract ive one, mediated 
by the pomeron exchange. 

The discussion in this paper has so far tacitly assumed a nucleon target. For nuclei, 
one extra effect to be taken into account is the longitudinal momentum transfer that is 
cut off by the nuclear form factor: 

do'-n-A-tqqA d&KN^qqN ,-,2 r 2 j-2 



dq\dz dq]_dz 



•FiKd, (4-5) 



where ( = q^/zzs, cf. ( |1.3| ), Fa\P\ = exp[— l/6R\t] and Ra is the mean square radius of 
the nucleus. We use Ra = 5.27 fm for platinum fl4| . In principle, one also has to take into 
account that the calculated initial momentum fraction distributions at the quark level are 
modified by nonperturbative hadronization corrections and effects due to the experimental 
acceptance. These effects lead to a certain broadening of the jet distributions compared 



the quark ones, as illustrated in Fig. 2 in Ref. [Q. For simplicity, we ignore them in what 



follows. We also ignore a tiny QED contribution of Coulomb scattering, calculated in |38 



The comparison of the calculated dijet momentum fraction distribution with the data 
14| in the transverse momentum range 1.5 < qj_ < 2.5 GeV is presented in Fig. [15|. The 
two solid curves correspond to the first two choices of the pion distribution amplitudes in 
Eq. ( f4.2p — asymptotic and "two-humped" , respectively. The dashed curve corresponds to 
the Chernyak-Zhitnitsky model evolved to the scale /i = 2 GeV. The overall normalization 
is arbitrary, but is the same for all three choices of the distribution amplitude. It is seen 
that experimental uncertainties do not allow for the separation between the distribution 
amplitudes <pi(z) and 03 (z) while the extreme choice <p2 {z) (see ( f4.2|) and Fig. jT2j) is not 



favored. This general conclusion is in agreement with the analysis in [[2"TH . 

Finally, one should not forget that the calculation involves a major uncertainty in 
the choice of the infrared cutoff parameter in the integration over the quark momentum 
fraction in the pion, Eq. (f4.1|). Our choice /iiR = 500 MeV corresponds to z' min = 1 — z' max = 
0.0625. The sensitivity to the IR cutoff is illustrated in Fig. [16] where in addition to our 
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Figure 15: The longitudinal momentum fraction distribution of the dijets with 1.5 < q± < 



2.5 GeV from the platinum target [14]. The two solid curves show the calculations with the two 
extreme pion distribution amplitudes in Eq. ( |4.2| ) — asymptotic and "two-humped" , respectively. 
The dashed curve corresponds to the Chernyak-Zhitnitsky model evolved to the scale \x = 1 GeV. 
The overall normalization (the same for all curves) is arbitrary. 
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Figure 16: Sensitivity of QCD predictions for the shape of the dijet momentum fraction dis- 
tribution to the choice of the infrared cutoff for the three different pion distribution amplitudes 
in Eq. (|4.2|) . The solid curves are calculated with z' min = 0.0625, the two dashed curves with 
z min = 0.0306 and z' min = 0.09, respectively. The normalization is adjusted to give the same 
integrated cross section in all cases. The data points are from [pi|]. 

"standard" choice (solid curves) we present calculations using /iiR = 350 MeV (z' min = 
0.0306) and /im = 600 MeV (z' min = 0.09), shown by dashed curves, for the three trial 
pion distribution amplitudes. The normalization of curves in Fig. |1^ is adjusted to give the 
same integrated cross section in all cases. We see that for the asymptotic pion distribution 
function 4>i(z) there is no effect at all (within the line thickness), the uncertainty is small 
for 03 (z), but more significant for 02 (z)- It is worth while to mention that the well-known 
strong scale dependence of the gluon distribution at small x (large energies) by itself leads 
to an effective suppression of the end-point contribution as compared to the factorizable 
part of the cross section. Indeed, the end-point contributions effectively correspond to large 
transverse distances and hence involve a (smaller) gluon distribution at low scales. This 
effect can be taken into account using the momentum fraction dependent factorization 
scale fi 2 F ~ q\(z' z') / (zz) , which may present an attractive alternative to the explicit 
cutoff. It is discussed in our letter [^] and yields results that are qualitatively similar to 
the calculation presented above. 

5 Conclusions 

We have presented a calculation of the leading contribution to the cross section of pion 
diffraction dissociation into two jets with large transverse momenta, originating from a 
hard gluon exchange. Our main result is that collinear factorization is violated. In techni- 
cal terms, the problem is due to pinching of singularities between soft gluon (and quark) 
interactions in the initial and in the final state. We have given a detailed analysis of 
this phenomenon using different techniques, and also explained how the structure of these 
singularities interferes with gauge prescriptions of the t-channel gluon propagators. Final 
expressions for the numerous contributions to the amplitude are collected in Appendix B, 
where in difference to the main text we use Ji's conventions |24| for the generalized parton 
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distributions. Using a realistic set of the generalized parton distributions [35] we find that 
the coherent dijet production is dominated by the gluon contribution. We also find that 
the real part of the amplitude is quite sizeable and cannot be neglected. 

It happens that the non-factorizable soft contribution to the dijet cross section imitates 
the shape of the asymptotic pion distribution amplitude. This implies that the dijet 
production is unlikely to yield fully quantitative constraints on the distribution amplitude. 
On the other hand, if one accepts that the pion distribution is close to its asymptotic form, 
as strongly suggested by the CLEO measurement [0 of the 77*7? transition form factor, 
then it appears that the QCD calculation explains the E791 data on the momentum 
fraction distributions of dijets very well. We find that the nonfactorizable contribution is 
suppressed compared to the leading contribution by a logarithm of energy so that for very 
large energies, in the double logarithmic approximation lng^ln£, collinear factorization 
is restored. This limit is achieved for energies about two orders of magnitude larger than 
the energy in the E791 experiment. 

The calculation presented in this work does not address specific nuclear effects, apart 
from the minor correction of the dijet distribution due to the nuclear form factor. For 
factorizable contributions of small transverse distances the dijet longitudinal momentum 
distribution is not expected to exhibit any significant A-dependence. The place where 
nuclear effects can play a role is by making the infrared cutoff introduced in ( f4.1| ) A- 
dependent. Indeed, it is plausible to assume that /ij R rises with the atomic number, e.g. 
~ A 1 / 3 , and hence nonfactorizable contributions become numerically suppressed due to 
color filtering of configurations with a large transverse size. The corresponding calculation 
goes beyond the scope of this paper. 

A short comparison with other approaches is in order. Our calculation is close in spirit 
to [20 , |21|j although the conclusions are different. A detailed explanation of the differences 
is given in the text. In short, the point is that we do not assume the light-cone dominance of 
the cross section from the beginning, but examine the light-cone limit carefully and argue 



that the approximation used in pO], [HJ breaks down for Glauber gluons (and quarks). In 
the double logarithmic approximation our result in ( 3.6(j| ) is similar to |T7J obtained using 
different methods. We, therefore, agree with the interpretation suggested in |TTJ that in 
the true diffraction limit, for very large energies, the dijet production can be considered as 
a probe of the hard component of the pomeron. We note, however, that this interpretation 
breaks down beyond the double logarithmic approximation and is not sufficient for the 
energy region of the E791 experiment. Finally, we have to mention an approach to coherent 
diffraction suggested in [18] that attributes hard dijets to a hard component of the pion 
wave function as in the original Brodsky-Lepage approach 0. This technique is interesting 
since it is most directly related to the classical view on diffraction B, but apparently 



complicated for the discussion of factorization. The general argumentation in [Qq] appears 



to be in contradiction with our explicit calculations. To our point of view, the coherent 



states formalism of Refs. 30 can be useful in this context. 
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Appendices 

A Generalized parton distributions in symmetric no- 
tation 

In the current literature two different sets of the generalized (off-forward, skewed etc.) 
parton distributions are used, as introduced in Refs. |23|, |24|, respectively. Our discussion 
in the text assumed the definitions Eqs. ( |2.1| ) and fl3.ip that correspond to the conventions 



of [ 23 1 . The transition to the notations of p4j] is straightforward and corresponds to the 



redefinition 

(1 + t)H q (y, = F q ( (x) - J^"(C - z)9(-l + C < x < C) , 
(1 + £)H g {y, = Ft(x) + ^(C - x)0(-l + C < x < C) • (A.l) 
The two sets of variables (x,Q and (y,£) are related according to 

,._x-C/2 C/2 

5 S 



a 1 - C/2 1 — C/2 

x -y + z C -JL (A2 ) 

y varies in the interval — 1 < y < 1. 

Rewriting of our results to the notations of |24| is straightforward. Using the crossing 
symmetry of the gluon coefficient function T gluon (x) = T g i uon (C — x), one can easily verify 
that 

i i 

'y + C 



.M g l U on ~ J dxJ 7 l(x)T g i nw (x) = \ J dy1-Lg(y,0 T glvLon 



(A.3) 



-1 

Note that the gluon distribution TC 9 (y,£) is a symmetric function, H g (y,£,) = 7~£g(—y,0- 
In the quark case we have 

i i 

quark ~ Jdx [j*{x)T q {x) + J*{x)T q {x)] = f dyH q {yii)T q fy^|j , (A.4) 

-1 

since quark and antiquark coefficient functions are related, T q (x) = —Tg(( — x), cf. (f2] 
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Figure 17: Gluon (left) and flavor-singlet quark (right) generalized parton distributions, for 
three different values of the asymmetry parameter £ = 0.008, £ = 0.012 and £ = 0.02 (from top 
to bottom) at the scale fj,p = 2 GeV, [|35|| , based on the LO MRST2001 parametrization 



In the numerical calculations in this paper we have used the set of leading-order parton 



distributions by Freund and McDermott based on the LO MRST2001 parametrization |36 
see Fig. [L7[ 



B Summary of Results 



Here we collect our final results. The hard coherent dijet production differential cross 
section is equal to 



da 



n-*qq 



dq\dtdz 



where we separate the quark and gluon contributions 

1 



<!J (/gluon ~\~ </quark) • 

7T 



(B.5) 



;b.6) 



The gluon contribution reads 
i i 

1 



-/gluon = - dyH g (y,0 / dz (j)«(z)l C F + J, 



, zz 
M — 7 + 1 



X 



1 / z z 



2N r \z> z' 



+ 



1 



y(z — z') — £(zz' + z'z) + ie y(z' — z) — i{zz' + z'z) + ie 



z' z'J 2N c z'z'J (y + £ - ie)(y - £ + ie) _ 

< 2 ™> { c *§ (I + ?) + 2^f (i + 0) + }■ < B - 7 > 
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The contribution in the last line of (|B.7[) corresponds to the pole contribution discussed 
in the text; the rest of the terms correspond to the naive light-cone limit £ — > £ — ze. The 
corresponding imaginary and real parts are 



-Im(J gk 

7T 



and 



zz 

+ ( — + 1 



zz 



1 / z z 

1 +^ - + 



2N r \z' z' 



y, r £,(zz'+z'z) 
9\ \z'—z\ 







\Z — Z\ 



" (Ml + J) + 2^) 



(B.8) 



1 1 



Re (Jgi U on) = dz' dy<f)n(z') lc F (J^ + ^ 



z z\ {2£H g (y,£) 







(y + 2 

2£(W 9 (y,£)-W s (£,£)-(y-£)F'(£,£)) 



(y - 5 



, ^z 
+ ( — + 1 



zz 



+ 1 + 



2; ^ 
2iV; \z~ l + l' 



x 2f (V?' 4- z'~) 11 
x z^zz + z z) _ _ , + ^ 



2£ (W ff (y,£)-W ff (£,£)) 



/ d^CO ^ C7f (J + p) (2£W;(£,£)ln 



l-£\ 2ft 9 (£,£) 



o 

+ — + 1 



£ 



1-^ 







(- 




\z' 





X 



In 



z — z 



|z' — z| — £(zz' + z'z) 



\z' — z| + £(zz' + z'z) 



1-e 



(B.9) 



where ^(£,£) = (d/dy)TC g (y,£,)\ y= £. Although the gluon distribution is generally not an 
analytic function at y — £, we have checked that the first derivative exists for the leading 
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order perturbative distribution ( |3.16| ) and one can prove that it also exists in the formal 
/ip — > oo limit. 

For the quark contribution we distinguish between the flavor-singlet contribution me- 
diated by the hard two-gluon exchange in the t-channel, and the separate contributions of 
u and d quarks: 

2£zz 

•/quark = ^ T (</ s ingl + Ju + J d) ■ (B.10) 

We also introduce the flavor-singlet quark distribution 



p=u,d,s 



(B.11) 



For the flavor-singlet contribution we obtain 



^singi = J dyH q (y,£) J dz' (j) n (z') C F + 1 



2y 



y 2 {z — z') 2 — t; 2 (zz' + z'z) 2 + ie ' 

(B.12) 



so that 



-Im ( J sing i) 

7T 



1 



X 



7-L 



£(zz' + z'z) 



Z — -2 



-£(zz' + z'z) 



\z — z\ 



(B.13) 



and 



} } / zz \ 2 y (w,(j/,o - ( 

R( '<-W : / dyy + 1) v y2(z _ z , )2 _\ 

1 



/ sign(y)S t (zz' +z' z) 



\z-z'\ 



+ J dz'^{z')C F (J| + l) In 




y\z - z') 2 - i\zz' + z'z) 2 
(z - z'f - i 2 {zz' + z'z) 2 



X 



Tin 



^(zz'+z'z) 



1 \ \z-z'\ 



^(zz' + z'z) 2 

i(zz'+z'z) 



\z'-z\ 



(z - z') 2 



(B.14) 



For the flavor-nonsinglet contributions we obtain 



1 1 

Ju = 2C F J dyH u (y,0 J dz'^(z') 
-1 

1 



z(2-z) 



+ 



zz'(y + £ — ie) z'{y — £ + ie 



+ 



z 



+ 



z — z 



A3 



2N C \zz' 2 (y + £-ie) z'z' 2 z(y(z' - z) - i{zz' + zz') + ie) 
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zz 



z(y(z' — z) — £(zz' + z'z) + ie) z'z' 2 (y(z — z') — £,(zz' + z'z) + ie) 

-(2m) (c/-^l + -L^) 5(y + 

V zz' 2J\ r zz' 1 ) 



(B.15) 



-lm(J u ) = 2C F J dz'(f) w (z' 
o 



v z zz 



+ 



Z — Z 



2N r \z'z' 2 z\z 



— z\ \ 1 — z — z J zz z 



1 „ ( Z(zx? + Jz) \ ^ / £(z^+zj) 



z — z 
z(2 - z) 1 z 



z z \z — z\ 



z — z' 



+2 CV v , y + 



zz' ' 2N r zz' 2 



w«(-£,0 



(B.16) 



and 



l i 

J d = 2C F J dy J dz'HdiyjO&Az') 

-1 o 

1 



zz'(y-£ + ie) z'(y + £-ie) 



+ 



z 



(z' - z) 3 



2N C \zz' 2 (y - £ -He) z'z' 2 z(y(z' - z) - £(zz' + zz') + ie) 
1 zz 



z(y(z' — z) — £(zz' + z'z) + ie) z'z' 2 (y(z — z') — £(zz' + z'z) + ie) 



z(l + z) 1 z 



+ (2ni) C F v ' 7 + 



zz' 2iV c zz 



/2 



«(y - 



(B.17) 



-lm(J d ) = 2C F [ dz'fa(z') \c F f^±^-H d ^,0 - ^7U(-U)) 

7T J I \ ZZ' Z' J 

£(zz' + zz' 



y - z) 3 



2N r V z'z' 2 z|z — z'\ V 1 — z — z' I 



\ ^J ^zz' + z'z) ^ 



z\z — z\ \ z' — z 
z(l + z) 1 z 



i(zz' + z'z 



zz 



z'z' 2 \z' — z 



H 



z — z 



-2 C 



+ 



zz' 2N r zz 



12 



(B.18) 



The expressions in the last lines of Eqs. (|B.15|) -( |B.18|) correspond to the pinch contribu- 
tions, see text. We do not present explicit expressions for the real parts of the u- and 
ci-quark contributions because the corresponding contributions to the dijet cross section 
are very small. 
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